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Abstract 

The charge transfer inefficiency (CTI) of the X-ray CCDs on board the Suzaku satellite (X-ray Imaging 
Spectrometers; XIS) has increased since the launch due to radiation damage, and the energy resolution 
has been degraded. To improve the CTI, we have applied a spaccd-row charge injection (SCI) technique to 
the XIS in orbit; by injecting charges into CCD rows periodically, the CTI is actively decreased. The CTI 
in the SCI mode depends on the distance between a signal charge and a preceding injected row, and the 
pulse height shows periodic positional variations. Using in-flight data of onboard calibration sources and of 
the strong iron line from the Perseus cluster of galaxies, we studied the variation in detail. We developed 
a new method to correct the variation. By applying the new method, the energy resolution (FWHM) at 
5.9 kcV at March 2008 is '^155 eV for the front-illuminated CCDs and '^175 eV for the back-illuminated 
CCD. 
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1. Introduction 

X-ray charge coupled devices (CCDs) have good spatial 
and energy resolution, and they have been the main de- 
tector for imaging spectroscopy in X-ray astronomy since 
the ASCA SIS (Burke et al. 1993). X-ray CCDs in or- 
bit, however, suffer from radiation damage. The dam- 
age causes the increase of the charge transfer inefficiency 
(CTI), which results in the degradation of the energy res- 
olution for two reasons: 1) the pulse height strongly de- 
pends on the position of an X-ray event, since the X-ray 
event loses more electric charges as the number of trans- 



fer increases, and 2) the loss of charge is a stochastic pro- 
cess, and thus there is a fluctuation in the amount of lost 
charge. In the case of the X-ray Imaging Spectrometer 
(XIS; Koyama et al. 2007) on board the Suzaku satel- 
lite (Mitsuda et al. 2007), the energy resolution at 5.9 keV 
was -140 eV (FWHM) at August 2005, and had degraded 
to ~200 eV at August 2006. 

The XIS is equipped with a charge injection (CI) struc- 
ture (Prigozhin et al. 2004; Bautz et al. 2004; LaMarr et 
al. 2004; Prigozhin et al. 2008) which lies adjacent to the 
top row of the imaging area. The CI structure allows us 
to inject a commandable amount of charge in a nearly 
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arbitrary spatial pattern. We can measure the CTI of 
each column precisely by the checker flag CI technique, 
and it is possible to correct the lost charge. The column- 
to-column CTI correction improves the energy resolu- 
tion greatly (Nakajima et al. 2008; Ozawa et al. 2009). 
However, we cannot correct the fluctuation in principle, 
and the degraded energy resolution cannot be fully re- 
stored even with the column-to-column correction. 

The CI can be used in another way to mitigate the 
effect of the radiation damage; the spaced-row charge in- 
jection (SCI) technique can reduce the CTI actively and 
improve the energy resolution. In the SCI technique, a 
charge is injected into CCD rows periodically. The in- 
jected charge fills the radiation-induced traps as a "sacrifi- 
cial charge" , and thus prevents some of the traps from cap- 
turing a signal charge produced by X-rays. Results based 
on ground experiments using the SCI technique with 
radiation-damaged CCDs have been reported (Tomida et 
al. 1997; Bautz et al. 2004), but no in-orbit experiment 
had been done. The Suzaku XIS operated with the SCI 
technique in orbit for the first time at August 2006, and 
the energy resolution was improved from ~200 eV to 
~140 eV at 5.9 keV (Bautz et al. 2007). The SCI has 
been a normal observation mode since October 2006. 

In this paper, we report the study of the CTI in the 
SCI mode, and present a new CTI correction method. 
Because the CTI depends on the distance between a pixel 
and a charge injected row, the calibration for the SCI 
mode becomes complicated. We developed a new method 
to correct the complex CTI of the SCI mode to improve 
the energy resolution further. We also found that the CTI 
increased with time even with the SCI, and our method 
can correct the time variation. All errors are at the Icr 
confidence level unless otherwise described. 

2. Spaced-row Charge Injection 

In the SCI mode of the Suzaku XIS, a charge is in- 
jected into every 54th row. The amounts of the injected 
charge into each pixel are equivalent to the X-ray energy 
of ^^6 keV and of '^2 kcV, for the front-illuminated (FI) 
and back-illuminated CCD (BI), respectively. 

The terminologies and notations for the CCD, CTI, and 
CI are the same as Ozawa et al. (2009). To make clear 
for further descriptions of the SCI, relevant terminologies 
and notations are summarized in table 1. 

The CTI in the SCI mode is assumed to consist of two 
components as CTI 1 and CTI 2 following Ozawa et al. 
(2009). The relation among readout pulse height {PH'), 
original pulse height generated by an X-ray of energy E 
{PHo at E), CTI 1 (ci), CTI 2 (cs), and transfer number i 
are formalized in equation 2 of Ozawa et al. (2009), where 
the detector coordinates are defined as ActX and ActY. 
The CCD consists of four segments (A, B, C, and D); each 
segment consists of 256 columns along ActX and has a 
dedicated read-out node. 

In the SCI mode, a radiation-induced trap is filled with 
a probability p when a sacrificial charge passes through 
the trap. The filled trap does not capture a signal charge 



produced by an X-ray. However, the electron once filling 
the trap will be re-emitted with the time scale r. Thus 
the probability of the trap holding the electron is propor- 
tional to p-exp(— f/r), where t is the time elapsed since 
the sacrificial charge passes the trap. Thus it is reasonable 
to assume that a pixel which is j rows away from its pre- 
ceding charge injected row has the ci value proportional 
to 1 — p • exp(— (5t • j/t), where St is time for one vertical 
transfer; 6t ■ j represents a time lag between the pixel and 
the charge injected row. We assume 6t ■ j /t is a small 
value, and hence I ~ p ■ exp(—St ■ j /t) is a linear function 
of j approximately. The charge is injected into every 54th 
row, and hence j = i mod 54 ^. We, consequently, can 
model the ci-i relation of the SCI mode as, 

Ci(*)=Cit + ^i^^-^-(*mod54). (1) 
54 

Equation 1 is a periodic sawtooth function as demon- 
strated in figure la. This ci-i relation generates the PH'-i 
relation in figure lb, and the shape is nicely reproduced by 
the ground experiments using the heavily damaged CCD 
as demonstrated in figure 15 of Tomida et al. (1997). 

The sawtooth distribution is also found in the XIS in- 
orbit. We measured the PH'-i relations using the onboard 
calibration source data in October 2006 with an effective 
exposure time of 1 Ms, and show the results in figures 2a 
and b. Figure 2 shows that the sawtooth of the BI sensor is 
shallower than that of the FI sensors. The results based on 
the data in February 2008 are also shown in figures 2c and 
d. Compared with October 2006, the sawtooth became 
deeper. The pulse height just after the charge injected 
rows changed more in the case of the BI sensor than the 
FI sensors. This difference might be resulted in part from 
the smaller amount of charge injected in the BI sensor. 

Our sawtooth CTI model can represent the complicated 
relation between PH' and i, and make it possible to con- 
vert PH'{i) to PHo with only three parameters, cit, cib, 
and C2. Our goal is to convert PH'{i) to PH^ by deciding 
the three parameters. 

3. CTI Measurement in Orbit 

3.1. Calibration Data 

To study the CTI in the SCI mode, we analyzed the 
data of the onboard calibration sources ^^Fe, the Perseus 
cluster of galaxies, and IE 0102.2—7219, whose properties 
are summarized in section 4 of Ozawa et al. (2009). 

XIS 2 suddenly showed an anomaly on November 9, 
2006, and it has not been operated since then. Although 
there is no direct evidence, the micro-metcoroid impact 
might have caused the anomaly ^. Thus only the data of 
XIS 0, 1 and 3 are studied. 

^ The modular arithmetic i mod 54 = j means that when i is di- 
vided by 54, it leaves j as the remainder, hence i — j is divisible 
by 54. For example, 111 mod 54 = 3. Strictly speaking, charges 
are injected in rows of j = 54 • n — 1 (n = 0, ■ ■ ■ , 18) and 1023 
for full window mode. Thus j should be (i + 1) mod 54. We 
nevertheless use this simplified description to avoid unnecessary 
confusion. We treated j properly for actual CTI measurement. 

^ See http:/ /www. astro. isas.jaxa.jp/suzaku/proposal/ao3/suzaku_td/. 
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Fig. 1. Our model of the "sawtooth" variation. We assume the ci-i relation show in (a), and it generates the PH'-i relation shown 
in (b). 
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(c) XIS 0(FI) Segment A Feb. 2008 
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(b) XIS 1(BI) Segment A Oct. 2006 
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(d) XIS 1(BI) Segment A Feb. 2008 
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Fig. 2. Pulse height of the Mm Ka line from the onboard calibration source as a function of i: October 2006 (a and b) and 
February 2008 (c and d). We show the results of the segment A in XIS 0, 1 as typical examples. X-ray events of grade 02346 were 
analyzed. Black and red marks are data before and after our new CTI correction, respectively. 
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Table 1. The list of notation. 

Notation Meaning 

i Transfer number. i=ActY + 1; ActY is a coordinate value where an incident X-ray generates a charge. 

PHo Original pulse height. 

PH'{0) Readout pulse height of a pixel at i — 0. It is equal to (1 — 1024 • C2) • PHq. 

PH'{i) Readout pulse height of a pixel at i. 

j Row number between a pixel and its preceding charge injected row. j = i mod 54. 

ci{i) CTI 1 of Ozawa et al. (2009) for a pixel at i. 

cit, cib The ci values at the peak and valley of the sawtooth (see figure la). 

C2 CTI 2 of Ozawa et al. (2009). 

—s{j) Slope for a pixel with j. It is equal to ci(i) • PHq. 

—St Slope for the tops of the sawtooth. It is equal to cn • PHo. 

— Sb Slope for the bottoms of the sawtooth. It is equal to Cib • PHq. 

P CTI depends on the pulse height as ci_2 oc PH~^ . See Ozawa et al. (2009). 



All data were acquired with the normal clocking mode 
and the full window option using the SCI. The editing 
mode was 3 x 3 or 5 x 5. Koyama et al. (2007) provide 
details of these modes. 

As mentioned in Koyama et al. (2007), a small fraction 
of the charge in a pixel is left behind (trailed) to the next 
pixel during the transfer. This charge-trail phenomenon 
changes the spatial extent of an X-ray event. All data 
were corrected for the phenomenon based on the in-orbit 
data ^, otherwise some X-ray events would be judged as 
grade 7, and the detection efficiency would decrease. 

We used the archival trend data of the calibration source 
obtained between August 2006 and March 2008, and the 
total effective exposure time is about 23.7 Ms. The obser- 
vations of the celestial objects are summarized in table 2. 

3.2. Determination of the CTI Parameters 

The procedures for the CTI determination are as fol- 
lows: 

• Step 1: deciding PH^ at 5.895 keV and 6.56 keV. 

• Step 2: measuring C2 for PH^ at 6.56 keV. 

• Step 3: measuring cu and cib for PHq at 5.895 keV 
or 6.56 keV. 

• Step 4- deciding the CTIs for any PH' values. 

In the case of the non-SCI mode, wc can measure the 
CTI of each column by the checker flag CI (Nakajima et 
al. 2008; Ozawa et al. 2009). Since the checker flag CI is 
a complicated operation, we have not used this technique 
in the SCI mode. We, therefore, measured the averaged 
CTI of each segment. 

In the normal analysis of the XIS data, both of single- 
pixel and multi-pixel events (grade and 2346 events; see 
Koyama et al. 2007) are used. To determine the CTI, 
however, we analyzed only the grade events; if we use 
the multi-pixel events, it is difficult to measure the CTI 
correctly, because the CTI depends on the amount of 
transferred charge, and the amounts of charge in each 
pixel comprising the multi-pixel event is different with 

^ The details about the charge-trail correction is shown in 
http: / / xmm2.esac.esa.int / external / xmm_sw_cal/icwg/ 
presentations / Suzaku_XIS .pdf 
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Fig. 3. PH' as a function of i obtained using the Fexxv 
Kcf line from the Perseus cluster. The result of the XIS 1 
segment C is shown as a typical example. 

each other. Wc will mention how to correct the data of 
grade 02346 events in step 4. 

Step 1: PHo at 5.895 keV and 6.56 keV 

We determined the PHoS of the FeXXV Ka line from 
the Perseus cluster (6.56 keV) by using the data of the 
first SCI observation on 29 August, 2006. Because it is 
impossible to measure the C2 value in August 2006, we as- 
sumed that value is equal to zero. Even if this assumption 
would not be reasonable, we can cover it by adjusting the 
PHo-E relation. Then the PH' ati^O (Pi?'(0)) is equal 
to PHo. 

Since the statistics are limited, we divided each segment 
to four regions along the ActY axis, and obtained the 
center pulse height of the iron line from each region. An 
example of PH' as a function of i is shown in figure 3. By 
fitting the data with a linear function of i, the PHo value 
was obtained as PH'{0) 

For the segments A and D, we obtained the PHo value 
of the Mnl Ka line (5.895 keV) using the data of the 
calibration sources in August, 2006. We obtained the 
PH'-i relation, but in this case, the i values are limited 
at around 900 because the calibration sources irradiate 
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Table 2. Log of Observations for Calibration. 



Obs. ID 


Observation time (UT) 


Exposure time 




Start End 


[ks] 


The Perseus cluster 


101012010 
101012020 
102011010 
102012010 


2006/08/29 18:55:07 2006/09/02 01:54:19 
2007/02/05 15:57:48 2007/02/06 14:30:14 
2007/08/15 12:40:49 2007/08/16 11:27:22 
2008/02/07 02:09:42 2008/02/08 10:30:19 


50.0 
43.9 
42.3 
41.6 


IE 0102.2-7219 



101005090 
101005110 
102001010 
102002010 
102003010 
102004010 
102005010 
102022010 
102006010 



2006/12/13 
2007/02/10 
2007/04/10 
2007/06/13 
2007/08/12 
2007/10/25 
2007/12/01 
2008/02/14 
2008/03/15 



18:53:16 
22:13:47 
10:35:08 
10:10:12 
05:21:09 
12:24:45 
19:25:40 
16:57:28 
05:43:27 



2006/12/14 
2007/02/11 
2007/04/10 
2007/06/14 
2007/08/13 
2007/10/26 
2007/12/02 
2008/02/16 
2008/03/15 



03:04:19 
19:30:14 
19:30:19 
03:31:19 
03:45:24 
09:00:14 
09:50:19 
03:10:24 
20:45:24 



28.2 
36.2 
18.1 
27.9 
39.5 
26.2 
24.8 
26.5 
28.2 



only the two far-end corners from the read-out node of 
the imaging area (Koyama et al. 2007). We extrapolated 
the relation to i = by using a linear function, the slope 
of which was fixed to a value smh calculated as follows. 
Because the slope represents an averaged — ci • PHo, the 
slope should be proportional to PH^~^. Then we calcu- 
lated SMn from the slope for the 6.56 keV data (spe) as 
SMn = SFc ■ (5.895/6.56)1-'^ ^ . (o.go)!"''. 

Step 2: C2 /or PH^ at 6.56 keV 

We obtained the PH'{0) values for all data of the 
Perseus cluster (table 2) with the same method of step 
1. The PH'-i relations are shown in figure 3. As is found 
in figure 3, PH'{Q) is decreased with time; PH'{0) be- 
comes different from PHq. We assumed the difference 
PHo — PH'{0) is attributable to C2. Assuming that C2 
is a linear function of time, the increasing rate of C2 for 
PHo (at 6.56 keV) is typically - 2 x 10"^ yr"! for the FI 
sensors and ~ 6 x 10^^ yr^^ for the BI sensor. 

Step 3: cit and cib for PHq at 5. 895 ke V and 6. 56 ke V 

For the segments A and D, we determined St and Sb 
for PHo at 5.895 keV by using the data of the cahbration 
sources. For the segments B and C, we measured those 
for PHo at 6.56 keV by using the Perseus data. We fitted 
the PH'-i relation with the sawtooth function shown in 
figure lb, and determined st and Sb. In the fitting, we 
fixed the PHo and C2 to the values obtained in step 1 
and step 2. Typical examples of the fitting arc shown in 
figures 4 and 5. Using st and Sb, we calculated cit and 

Clb- 

In the case of the FI sensors, the cn and cib values for 
PHo at 5.895 keV at August 2006 are typically < 10"'^ 
and '--^ 3 X 10"^, respectively. The increasing rates are 
^ 1 X 10^^ yr^^ and ~ 5 x 10^® yr^^, respectively. In 
the case of the BI sensor, the parameters cn and cib at 
August 2006 are - 6 x 10"^ and 1 x 10"^ , and the 
increasing rates are ~ 5 x 10^® yr"^ and 8 x 10"^ yr"^. 
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Fig. 4. Best-fitting result of the sawtooth function for the 
Mm Ka line of the calibration source. The result of the XIS 
segment A in February 2008 is shown as a typical example. 
Only grade events were used. 

respectively. Typical examples of the time evolution of 
the Ci values for PHo at 5.895 keV are shown in figure 6. 

The CTI of the BI sensor is larger, and increases more 
rapidly than that of the FI sensors. It might be the reason 
that the amount of the injected charge to the BI sensor is 
less than that of FIs. 

Step 4- CTI for any PH' values 

For an event of PH', we calculated 
the CTI with an equation of ci,2 x 
{PH'/{PHo at 5.895 kcV or 6.56 keV)}-'^. Here, 
we assumed the same PHo-ci,2 relation as that of the 
non-SCI mode and used the same (3 (typically ~ 0.25; 
Ozawa et al. 2009). 

Since we obtained the CTI for all pulse-height values, we 
can now correct the CTI of multi-pixel events (grade 2346 
events). By using the Mnl Ka line of the calibration 
sources or the FeXXV Ka line from the Perseus clus- 
ter, we examined the CTI correction for the grade 02346 
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Fig. 5. Best-fitting result of the sawtooth function for the 
Fexxv Ka line from the Perseus cluster. The result of the 
XIS segment C in February 2008 is shown as a typical ex- 
ample. The white and black circles show the pulse height of 
regions with smaller and larger j values, respectively. Only 
grade events were used. 
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Fig. 6. Time history of cit and cib for PHo at 5.895 keV. 
The segment A of XIS and 1 are shown as typical examples. 

events. Figure 7 shows the Hne centroid of the grade 
and grade 02346 events as a function of time. While the 
line centroid of the grade events is temporally constant, 
that of the grade 02346 events increases with time. This 
means that the CTI parameters for the grade 2346 events 
are smaller than those for the grade events. The CTI 
of the grade 2346 events can be smaller, because the split 
charge acts as the sacrificial charge. 

We then fine-tuned the parameters cn, cib, and C2 by 
multiplying a common time-independent factor. We de- 
termined the factor of each segment so that the line center 
of the grade 02346 events becomes temporally constant. 
The factors are typically 0.9. Figure 7 also shows that 
the result of the grade 02346 data corrected with the fine- 
tuned CTI parameters. The adjusted pulse height be- 
comes constant. However, as a result, the value becomes 
different from PHo measured in step 1. We regarded this 
adjusted pulse height as PHo at 5.895 keV, and deter- 
mined the PHo-E relation. 

To check the CTI correction in the low-energy band, we 



Fig. 7. Time history of the line centroid of Mn I Ka line. We 
show the result of the XIS 1 segment A as a typical example. 
The triangle and circle marks show the data of grade and 
02346 after the CTI correction with the parameters based 
on only the grade events. The rectangle marks show the 
data of grade 02346 events after the CTI correction with the 
fine-tuned parameters as described in the text. 
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Fig. 8. Time history of the line centroid of Ovill Ka 
(0.653 keV). We show the result of the XIS 1 segment C as 
a typical example. The results are based on the grade 02346 
events corrected with the fine-tuned CTI parameters as de- 
scribed in the text. 

applied the fine-tuned CTI parameters to the grade 02346 
data of IE 0102.2-7219. We fitted the spectrum with the 
empirical calibration model^. Figure 8 shows the pulse 
height of the O VIII Ka line after the CTI correction. We 
can see that the line centroid is constant, which supports 
the validity of our method. 

4. Energy Scale Uniformity and Resolution in 
the SCI Mode 

The pulse height of the Mn I Ka line after our new CTI 
correction is shown in figure 2. The sawtooth structure 
disappeared, and hence our new method greatly reduces 
the variation of the pulse height. Comparing between 
October 2006 and February 2008 shows that the CTI vari- 
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(a) 



ation with time is also corrected properly. 

The PHo-E relation for each segment is assumed to be 
the same function form with a similar fine-tuning pro- 
cess of the absolute energy as Ozawa et al. (2009): a 
model of two slopes crossing at the energy of the Si-K 
edge (1.839 keV) with the same ratio of the two slopes as 
that obtained in the ground experiments (Koyama et al. 
2007). 

The time histories of the energy resolution of the cor- 
rected data are shown in figure 9. At the high energy 
(5.895 keV), the energy resolution of XIS at March 2008 
is improved from ~160 eV to ^^155 eV by the sawtooth 
correction. On the other hand, the energy resolution of 
XIS 1 at March 2008 is not improved in spite of the saw- 
tooth correction, and stays at ~175 eV. At the low energy 
(0.653 keV), no clear effect of the sawtooth correction is 
seen in both FI and BI. The energy resolutions at March 
2008 are ~53 eV (XIS 0) and ~62 eV (XIS 1) indepen- 
dently of the correction. 

As we mentioned in section 2, the sawtooth of the BI 
sensor is shallower than that of the FI sensors. We think 
it causes the difference of the effect of the sawtooth cor- 
rection between BI and FI. We speculate, at the low en- 
ergy, the effect of the sawtooth correction was not seen 
because the original degradation of the energy resolution 
was small. In fact, the energy resolutions at O VIII Ka are 
almost same between the SCI and non-SCI mode. 

The results of our new method to correct the sawtooth 
structure has been implemented to the software package 
released by HEASARC^ since HEAsoft version 6.3. Now, 
all of the XIS data of the SCI mode after the processing 
version 2.0 are corrected by the sawtooth method. 
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Fig. 9. Time history of the energy resolution (FWHM) at 
Mm Ko (5.895 keV) (a) and at Oviii Ko (0.653 keV) (b). 
The averaged value of the segments A and D are shown for 
each sensor. The results are based on the grade 02346 events 
corrected with the fine-tuned CTI parameters as described in 
the text. 
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